Flat lens without optical axis: Theory of imaging 
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We derive a general theory for imaging by a fiat lens without optical axis. We show that the 
condition for imaging requires a material having elliptic dispersion relations with negative group 
refraction. This medium is characterized by two intrinsic parameters a and k. Imaging can be 
achieved with both negative and positive wave vector refraction if a is a positive constant. The 
Veselago-Pendry lens is a special case with a = 1 and k = 0. A general law of refraction for 
anisotropic media is revealed. Realizations of the imaging conditions using anisotropic media and 
inhomogeneous media, particularly photonic crystals, are discussed. Numerical examples of imaging 
and requirements for sub-wavelength imaging are also presented. 
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Since antiquity, the positive index of refraction of con- 
ventional materials has required the use of curved sur- 
faces to focus light. However there has been a continuing 
quest for lenses with flat surfaces as they confer a variety 
of advantages. Notable examples are the Fresnel lens and 
Maxwell's fish eye lens yj, |2| ■ The Fresnel flat lens uses a 
gradient-index n{x) material (GRIN) 0,0,0 an d hence 
must possess an optical axis, i.e. is not translationally 
invariant along the surface (x-axis). 

The concept of negative refraction has led to new fun- 
damental approaches as well as applications in optics 
0. In 1968 Veselago pointed out that for a mate- 
rial with refractive index n = — 1 a flat surface would 
focus light. A decade later, Silin discussed a more 
general case and obtained the lens equation of a flat slab 
with negative elliptic dispersion. Pendry's recent anal- 
ysis demonstrating the possibility of sub-wavelength 
resolution with a flat slab of such materials, as well as the 
experimental realization of the so-called left-handed ma- 
terials usin g co mposite media E3 an d photonic crystals 
(PhCs) mi |12| led to renewed interest in the unique elec- 
tromagnetic properties of these artificial materials. The 
unique property of these flat lenses is the lack of optical 
axis. This type of flat lenses can be realized in a PhC us- 
ing negative refraction [T^ | and has been demonstrated 
in microwave experiment |l4j| . Super-resolution imag- 
ing through single negative index media are also studied 
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There is however no complete theory 
of imaging by a flat lens which properly describes the 
various features observed in the experiments and in nu- 
merical simulations. 

In this paper, we present a general theory of imaging 
by a flat lens without optical axis, resulting in a pro- 
posal for a new material which is translationally invari- 
ant along the surface and has anisotropic refractive index 
n(9). Defining an optical phase condition for imaging, we 
show that the condition requires a material having Ellip- 
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FIG. 1: Wave vector refraction for (a) a single interface be- 
tween the vacuum and an ENGRM and (b) an ENGRM slab. 



tic dispersion relations with Negative Group Refraction 
(ENGRM) at the operating frequency. The ENGRM is 
defined by two parameters: the anisotropy parameter a 
which is a measure of the ellipse eccentricity, and a phase 
factor k which determines the center of the ellipse. Re- 
fraction laws for the wave vector and the group velocity 
for anisotropic media are derived. The theory shows that 
"perfect" images can be obtained, and imaging is possi- 
ble with negative as well positive refractive indices. The 
required conditions and consequences of flat perfect and 
imperfect lenses made of ENGRMs are derived from a 
generalized Fermat's principle. Realization of flat lenses 
using PhCs is discussed. In real PhCs, a is itself angle- 
dependent in most cases, leading to some limitations for 
image formation. Many of the proposals for PhC lenses 
are contained in the present theory. The theory is also 
applied to recent experiments using PhCs and is shown to 
successfully describe some intriguing aspects of the data. 
Numerical simulations are presented that describe visual 
details of the image formation in PhCs. 

We first consider the case of a single interface which is 
along the x-axis and at z = as shown in Fig. QJa). The 
wave vector k of an incident plane wave from vacuum 
towards the interface makes an angle 6 with the z-axis. 
Here we consider the case that for any incident k, there 
is only one refractive wave vector k r . Continuity along 
the interface requires that k rx = k x . A point source is 
at x — and z — — u with u > 0. We now consider 
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another point in the ENGRM at x = Ax and z = v with 
v > 0. The phase difference between these two points 
for each incident k is $ = $ 2 + "^with <E> Z = k z u+ 
k rz v and <& x = k x Ax. According to Fermat's principle, 
the formation of an image would require that the total 
phase <E> be stationary. This restriction can be relaxed 
for a flat interface. Since the surface is flat, any point 
of a finite-size source will be imaged to a point with the 
same x-coordinate. Thus Ax = and Q x = 0. One only 
needs to consider <£> z . A generalized Fermat's principle 
states that an image will be formed if the phase $ z is 
stationary, d$ 2 /dfc z = 0. The lens equation for a single 
interface is 



with a material constant 



-dk rz /dk z 



(1) 



(2) 



This constant a should be positive and independent of 
the incident angle for a focus without aberration. Thus 
the following rule for the wave vector refraction at the 
interface must be obeyed 



k — 



(3) 



Here k is the integration constant which is also an in- 
trinsic property of the ENGRM. Since in the vacuum 
k\ + k 2 x = fc 2 = Wg/c 2 , the equi- frequency surface (EFS) 
of the medium is 



'■{n-k rz ) 2 +k 2 rx =LU%/C 2 



(4) 



Note that this elliptic dispersion exists only at the op- 
erating frequency loo [Tflj . In the neighborhood of this 
frequency uj ~ u; , (k — k rz ) 2 /n 2 + k 2 x /n 2 = uj 2 /c 2 with 
a = n z /n x ~ a. The deduction of Eq. Q from an 
elliptic EFS requires negative group refraction 



(k r — KZ,) ■ Vk r LU < 



(5) 



together with the causality condition duj/dk rz > that 
the energy must flow away from the interface which 
is illustrated in Fig. [21 The ray vector [2(j which 
represents the direction of group velocity in a medium 
with elliptic EFS is s = (— ersinax + cosaz)/n x a while 
k r — kz =n x (sinax — trcosaz) with sin a = n" 1 sin# and 
the incident angle in the vacuum. 

The Snell's law for wave vector loses its meaning in 
this anisotropic medium. However there is a law for 
group refraction, which is tan/3 = — (rtana = —a(n 2 — 
sin 2 sin 9. At the operating frequency ujq of the 

flat lens, n x = 1, u{ujq) — a, the group refraction law is 
simply 



tan j3 



-crtanS. 



(6) 



Thus —a can be regarded as the effective refractive index. 
Note that this refraction law is very general and is valid 




FIG. 2: Group velocity for an ENGRM with elliptic EFS and 
(k r — kz) ■ Vk r ^ < and (a) negative or (b) positive wave 
vector refraction. The dashed arrow is k r — kz. (c) Ray 
diagram of the group velocity for imaging by an ENGRM flat 
lens with u + v = ad. 



for any a as defined by Eq. J2J ■ The proof is the follow- 
NS, tan/3 = = = -vtane. Here 

we used k rx — k x . 

The space inside the medium can be considered to 
be "optically stretched" in the z-direction by a factor 
of a . The waves inside the medium are represented 
by tP(y') — exp(ik x x — ik z z' + in' z') with z' = az and 
k' = k/ct. Thus in the case that k^O, the waves at this 
frequency will experience a periodicity in the ^-direction 
with length scale 27t/|k| since the waves are modulated 
by exp(iKz). This periodicity is the effective periodicity 
and may not be the actual physical periodicity such as 
the lattice spacing of a PhC. 

We now consider the imaging of an ENGRM slab with 
thickness d as shown in Fig. db). The first surface is 
at the origin and the second at z = d. A point source 
is placed at z = —u. We consider another point outside 
the ENGRM slab at z = d + v. If the value of v satisfies 
the following lens equation 



u + v = ad, 



(7) 



the phase difference <& z = k z (u + v) + k rz d = nd is sta- 
tionary and an image will be formed without aberration. 
Note that for an ordinary thin curved lens with focal 



length /, the lens equation is u 



r 



With 



given value of a and slab thickness d, the distance be- 
tween the object and the image outside the slab is fixed, 
u + v + d = (1 + a)d. Once the location of the object is 
fixed, the position of the image is also fixed, no matter 
where the slab is placed! 

We remark that a flat slab can focus with negative as 
well as with positive wave vector refraction. The neces- 
sary condition for flat lens imaging is negative group re- 
fraction. In general, there is not much meaning of Snell's 
law for wave vectors for anisotropic media. Only the re- 
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fraction law for the direction of group velocity Eq. © 
is meaningful and does correspond to ray diagrams. The 
Veselago-Pendry lens is a special case of our flat lens with 
a = 1, K = and is the only case that the ray diagram 
is applicable for the wave vectors. 

The lensmaker's formula Eq. (JIJ or only provides a 
necessary condition for imaging. An additional constraint 
is required for the formation of a "perfect" image, viz. 
the surface reflection coefficient r should vanish, r = 0, 
for far field and should diverge, r — > oo, for near field. 
These set the conditions for a perfect flat lens without 
optical axis. The requirement of no reflection r = is 
not essential for a slab to focus far field. The principal 
conditions for the flat lens to focus light are that Eq. ffil 
and \r\ <C 1 should be satisfied for all or at least a large 
range of incident angles. The presence of reflection will 
make the image dim and may give rise to multiple images. 
The indefinite indices medium |2l| can have an elliptic 
dispersion relation, but in general will not satisfy the flat 
lens imaging conditions. 

The field inside the ENGRM flat lens can also be de- 
scribed by a partial differential equation [a~ 2 (iK + d z ) 2 + 
d 2 + (ui/c) 2 ]^ — 0. The medium can be regarded as a 
distorted space with the metric ds 2 = dec 2 +a 2 dz 2 . Thus 
there is an analogy with gravitational lensing where re- 
fraction occurs due to the warping of space caused by 
general relativistic effects, such as in the vicinity of a 
massive object. It is worth noting that general relativis- 
tic effects are observable due to optical refraction. 

In the rest of the paper, we will focus on S'-polarized 
electromagnetic waves with the electric field in the y- 
direction. In general, the magnetic permeability /i of the 
ENGRM should be a tensor, and its relationship to the 
other material parameters a and k can be deduced. Here 
we assume that effective indices e and /i can be used to de- 
scribe the medium of interest. Due to its symmetry, the 
ellipsoid can always be transformed to its principal axes 
2J. For the S-polarized waves, the H field is in the xz- 
plane, only [i x and \i z will enter our discussion. Further- 
more, fj, x plays a more active role than /i z since the sur- 
face reflection coefficient is r = (^ x k z — k rz )/ ([i x k z + k rz ). 
Thus to achieve perfect imaging for far field, the effective 
permeability of ENGRM should be 



fj, x = n/ko cos ( 



(8) 



Only in the case k = that \x x = — a is isotropic while for 
k ^ 0, fj, x has to be anisotropic and diverges at k x = k$. 
For k ^ 0, any finite \i x at k x = kg results in reduced 
transmission. Notice that in PhCs including metamate- 
rials, the effective (j. x should be anisotropic in general. 

To see how the propagating waves are focused, we con- 
sider a S'-polarized point source — (i/4)iJQ 1 ^(/c |r + uz\) 
which is centered at z — —u while the two surfaces of 
the ENGRM slab are at z — and z — d. If \x x takes 
the form given by Eq. JSJ), there will be no reflection 
for propagating waves and the transmission coefficient 
through the slab is T = exp(ik rz d). The transmitted far 




FIG. 3: Field intensity of two sources ( Jo(fcor)) imaged by an 
ENGRM flat lens with a = 0.5, k = 0, and thickness d = 10. 
The two white lines indicate the surfaces of the lens. We set 
fi x = —0.5 for the ENGRM for perfect transmission. 



field is 



£4 r (r) = -ie lKd f f^ e tt.[*+»-(i+')<fl cos(k x x). (9) 
Jo 2n 

Note that k z = ko cos 9 and k x = ko sin 9. The image 
outside the slab is located at z — (1 + a)d — u. One can 
see that except for a global phase exp(iKd), the far field 
of the image is exactly that of the source. The far field 
inside the slab is 



-te 



*' 2 d9 



-<M«-«0 cos^a) (io) 



for < z < d. Another perfect image is formed inside the 
slab at z = u/a. The first surface acts as a mirror with 
an extra phase exp(iKz). The space inside the ENGRM 
slab is optically stretched by a factor of cr _1 . A flat lens 
imaging is shown in Fig. |3J 

Evanescent waves and sub-wavelength imaging For 
evanescent waves k x > ko, the complex extrapolation of 
Eq. (gj gives 



iaq 



(11) 



with q = (k 2 — fcg) 1 / 2 . The transmission coefficient 
through the flat lens is T = (1 — r 2 ) exp(ik rz d)/[l — 
r 2 exp(2ik rz d)] with r = [{[i x — <r)q + iK\/[(^, x + a)q — in]. 
In order to amplify evanescent waves, a singularity must 
exist in T to compensate to certain extent the decay 
of evanescent waves in the vacuum exp[ — q(u + v)]. In 
the case of single- interface resonance [3, |24|, r — * oo, 
Texp[— q(u + v)] = exp(iKd), all the evanescent waves 
will be amplified only if 



in/q. 



(12) 



The images both inside and outside have the same sub- 
wavelength features of the source. If fi x couldn't take 
the above form, sub-wavelength imaging is still possi- 
ble. At the so-called overall resonance condition [24| . 
1 — r 2 ex ] p{2ik rz d) = 0, thus r = ± exp(— ik rz d), one gets 



/'.< 



(-a + iK/q)t&nh Tl (o-q- in)d/2. (13) 
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For k 7^ 0, fj. x and fx x are all complex with negative real 
parts. One notices that H x {q) is flatter than fJ,~{q) as 
functions of q. In terms of effective indices, the existence 
of surface modes u>o(k x ) requires that JJ/iz be negative for 
5-polarization. To amplify evanescent waves, the u>o(k x ) 
curve must be very flat |2 1| . This is equivalent to say that 
the curve ri x (q) should be flat so that the contribution 
of evanescent waves are constructive for certain window 
of q. The closer fx x is to a4~(<z), the more sub- wavelength 
features the images will have. Note that for large nd, 
the evanescent waves are located in the vicinity of the 
interfaces for constant fi x . 

Realization in real materials We note that no known 
natural material is found to have negative dispersion. 
Negative group refraction can be achieved in periodic or 
quasi-periodic media with nonzero k 22]. For most me- 
dia, Eq. @) and ijjjjj) are unlikely to be fulfilled, thus a 
perfect flat lens is unattainable. Realization using real 
materials is discussed next. 

Ordinary material may not have the EFS described by 
Eq. (0J for all k x . For an uniaxial crystal, since the group 
refraction is positive, thus k rz — n — ak z with n = and 
a < 0. There will be no focus though negative refrac- 
tion for certain incident angles can be easily achieved 
|23j . Instead a virtual image will be formed satisfying 
the same lens equation. However, a PhC could have both 
an elliptic EFS and negative group refraction for certain 
windows of k x at certain frequency. Here we use some 
general features of PhCs to explain the imaging mecha- 
nism. To this end, we consider the first band of a square 
lattice PhC of unit lattice spacing. For small k x along 
the YM direction, the dispersion can be approximated 
elliptically as k rz ~ k — aok z with appropriate constants 
k and do (see the dashed line in Fig. At ujq = uj u , 
the upper limit of all-angle negative refraction (AANR) 
in the first band k = v^tt — fco(l — 00) ■ Since for 
AANR with ujq away from oj u , do will be smaller, we 
only concentrate here on the discussion for ujq = lu u . 
To enhance the transmission, the slab thickness can be 
selected to satisfy the Fabry-Perot resonance condition 
(k — gqlo 1 1 c)d = (mod 2tt). The slab is on resonance 
only for k x — > 0. For large k x , k rz ~ \/2iv — <r 1 ^ 2 k z 
due to the C^ v symmetry of the EFS. A desired large 
transmission would require /j, x to diverge. Thus the im- 
age location is bounded a < [u + v)/d < o~ and the 
Fourier components of the object with k x — > ko will be 
partially lost. This leads to a so-called self-collimation 
effect for small a®. Far field images both inside and out- 
side the slab with reasonably good quality will be formed. 
However, the image inside the PhC is stretched and mod- 
ulated due to the partial transmission besides the Bloch 
wave modulation. Even with the removal of Bloch wave 
modulation, the image inside a PhC will be visible only 
for large thickness. For evanescent waves at ui = co u with 
k x —* ko, one has = — cr ^ 2 and — > 00 while for 
k x S> ko, /Lt^(fe x ) ~ fj, x (k x ) — > — (Jq 1 - To have substantial 
sub-wavelength feature of the image, the real part of [i x 
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FIG. 4: (a) EFS of the TE modes of a square lattice PhC 
calculated using plane wave expansion, (b) Far-field H z of 
a point source across an ENGRM slab with dispersion rela- 
tion krz = K — ak z with k = 3.2465 and a = 0.08 which 
approximates the actual EFS of the PhC (dashed line in the 
upper panel). The permittivity e x = 2.2 is used. Note the 
modulated field inside the homogeneous ENGRM. Due to 
impedance mismatch, the images are not "perfect" unlike Fig. 
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must take large negative value for small cro- An example 
is shown in Fig. which is very similar to Fig. 5 in Ref. 
[T^j l . The image is confined in the vicinity of the second 
surface of the PhC slab due to a small value of cro- 

In this paper, we have discussed the group refraction 
of an anisotropic medium characterized by two materi- 
als parameters a and k. When a is a positive constant, 
the anisotropic medium is an ENGRM and a flat slab of 
such medium can be used as a focal lens without optical 
axis and leading to images free of aberration. Our theory 
of flat lens imaging is a generalization of the Vcsalago- 
Pendry perfect lens and beyond Silin's formula. The 
theory is valid for both real and virtual images. This 
flat lens theory leads to a clear understanding of neg- 
ative group refraction and flat lens imaging in electro- 
magnetism 0, ^J, acoustics [25|, and electron waves. 
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Numerical simulations of homogeneous ENGRMs and of 
PhCs were carried out supporting our theory. 

The lack of optical axis for the flat lens has very broad 
applications and confers important advantages in optics. 
Clearly there are no aperture restrictions. There are 
two characteristics that might be viewed as limitations. 
The lensmaker's formula clearly dictates upper bounds to 
u,v < ad for real image formation. The other limitation 
is that magnification is always unity. For any medium the 



lens properties a, n, and also the working frequency ujq 
can be obtained by inspection of the EFS. To have high 
quality image, a should be constant and close to unity. 
The present theory can be used to design tailor-made 
flat lenses. Extension of our theory to three-dimensions 
is straightforward. 

Work supported by the National Science Foundation 
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